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ABSTRACT 



Context. Accreting X-ray pulsars provide us with laboratories for the study of extreme gravitational and magnetic fields, hence 
accurate descriptions of their observational properties contribute to our understanding of this group of objects. 

Aims. We aim to detect a cyclotron resonance scattering feature in the Be/X-ray binary Cep X-4 and to investigate pulse profile and 
spectral changes through the outburst. 

Methods. Spectral fitting and timing analysis are employed to probe the properties of Cep X-4 during an outburst in 2002 June. 
Results. A previously announced cyclotron feature at 30.7 keV is confirmed, while the source shows spectral behaviour and luminosity 
related changes similar to those observed in previous outbursts. The long-term X-ray lightcurve shows a periodicity at 20.85 d, which 
could be attributed to the orbit in this Be system. 

Key words. X-rays: stars -X-ray binaries: cyclotron lines -stars: magnetic fields -stars: pulsars: individual: Cep X-4, GS 2138-1-56 
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1. Introduction 

Cep X-4 was discover ed as a transient source in 1972 
June/July with OSO 7 dUlmeret all Il973l) . In 1988 March 
it wa s detected with Ginga during a month long X-ray out- 
burst (' Makino & Ginga Teamll988ah and labelled GS 2138-H56. 
During this outburst, where the source reached a maximum in- 
tensity of lOOmCrab, pulsations at 66.2490 ± 0.0001 s were 
detected (iKoyama et alJ ll99U iMakino & Ginga TeamI [l988bl) 
as well as a cyclo tron resonance feature at 30.5 + 0.4 keV 
( lMiharaetalJll99lh . leading to its classification as an X-ray bi- 
nary pul sar. Using the ROSAT Position Sensitive Proportional 
Counter, ISchulz et al.l ([19951) detected Cep X-4 in quiescence 
in 1993 January and in outburst in 1993 June. During out- 
burst, the pulse period was measured at 66.2552 + 0.0007 s. 
This observation lead to a refinement of the source position, 
and a Be star at azooo = 2lh39"30:6 and ^ 2000 = -H56°59a2r 9 
was proposed as the optical counterpart dRoche et al.l Il997h . 
Long slit optical spectroscopy of this proposed counterpart 
jBonnet-Bidaud & Mouchetl Il998l) confirmed Cep X-4 as a 
Be/X-ray binary with a distance of 3.8+0.6 kpc. In 1997 July and 
August, another outburst took place, this time observed by the 
Burst and Transient Source Exp eriment (BATSE) an d the Rossi 
X-ray Timing Explorer (RXTE). IWilson et alJ (11999 ) measured 
the pulsar spin down rate and placed constraints on the orbital 
period using BATSE and RXTE data from the 1993 and 1997 
outbursts. The most recent outburst, occuiiing in 2002 June and 
lasting approximately one month, was observed by RXTE only. 
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In this paper, we describe the spectral and timing analysis 
of pointed RXTE observations during the 2002 June outburst. In 
Sect. 2 we present the observations and describe the reduction 
techniques. Section 3 presents a possible orbital period in the 
long term lightcurve while in Sect. 4 we fit X-ray spectra from 
the 2002 outburst. In Sect. 5 the pulse profiles are introduced. 
Section 6 presents a discussion and interpretation of the data, 
while the results are summarised in Sect. 7. 

2. Observations and Data Reduction 

Figure □ shows the RXTE All Sky Monitor (ASM) X-ray 
lightcurve of the month long outburst of Cep X-4 during the sec- 
ond quarter of 2002. Observations with RXTE's pointed instru- 
ments were performed during the outburst (aiTows in Fig.[T]and 
Table[B. 

To increase the signal-to-noise we extracted data from 
only the top layer of the Proportional Counter Array (PCA 
iJahoda et al.l 120061) . As most of the source photons, but only 
half the instrumental background, is detected in this layer it is a 
good choice for a relatively weak source such as Cep X-4. These 
data are Standard 2 mode data with 16 s time resolution and 128 
channel energy resolution, and were employed to generate phase 
averaged spectra in the range 3.5-20 keV. We added systematic 
errors of 0. 1 % of the count rate in quadrature to all spectral bins 
across the energy range, as indicated by near contemporaneous 
fits to the Crab spectrum. For pulse profiles and timing analysis, 
top layer Good Xenon data with a time resolution of 0.03 s were 
used. 

For all spectra we used the E_8jUj'_256JD XlF mode High 
Ener gy X-ray Timing Experiment (HEXTE iRothschild et al.l 
Il998l) data, which has a temporal resolution of 8 //s. In all cases. 
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Table 1. KXTE Observations of Cep X-4 during 2002 June out- 
burst. Luminosities are in the 2-10 keV range and are calculated 
using a distance of 3.8 + 0.6 kpc dBonnet-Bidaud & Mouche^ 

[THi. 



Obs ID 


MJD 


Exposure 


Luminosity 


Period 


70068-11- 






lO^^^ergs"' 


s 


01 -UO 




oOO 


1.40 


00.41(0.05; 


01-01 


52442.1 


816 


1.06 


66.36(0.04) 


01-02 


52444.7 


880 


0.87 


66.25(0.04) 


02-01 


52450.1 


1552 


0.41 


66.35(0.02) 


02-000 


52450.7 


19600 


0.35 


66.30(0.01) 


02-00 


52451.0 


1456 


0.31 


66.33(0.03) 


02-02 


52451.3 


8432 


0.32 


66.30(0.01) 


02-03 


52451.7 


7088 


0.27 


66.296(0.004) 
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Fig. 1. The ASM lightcurve of Cep X-4 during the 2002 June 
outburst, re-binned to a resolution of 1 d and taking only one-day 
averages comprised of 10 or more dwells into account. Arrows 
indicate the dates of pointed RXTE observations. 

the signal-to-noise ratio was increased by adding together data 
from both HEXTE clusters. Data above 70keV were not used 
due to rapid deterioration of the signal-to-noise at higher ener- 
gies. 

Analysis was done with HEADAS version 6.0.4 and spectral 
fitting with XSPEC version 1 1.3.1 w (,Ai-naud,1996,) . 

3. Long-term lightcurve 

The archival RXTE-ASM data from MJD 50090 to MJD 53873 
were filtered on the criterion that one-day averages should con- 
sist of 20 or more dwells to make a long-term lightcurve of time 
resolution one day. Without losing too many data points, this fil- 
tering gets rid of the very noisiest points, but includes data from 
both the 1997 a nd 2002 outbursts of Cep X- 4. Figure|2]shows the 
Lomb- Scargle (lLomblll976t IScargldl 19821) periodo gram of this 
X-ray lightcurve. To determine the significance levels of peaks 
in the periodogram, 10000 Monte Carlo white noise simulations 
were generated with the same mean, variance and sampling as 
that of the ASM data. The 90% and 99% significance levels are 
shown in Fig. |2] A peak corresponding to a period at 20.85 d, 
with an error of 0.05 d is found above the 90% significance level. 

Apart from this periodicity and red noise at low frequen- 
cies, no other significant long-term periodicities are found in 
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Fig. 2. Lomb-Scargle periodogram of the ASM lightcurve of 
Cep X-4 

10-lOOOd interval. This periodicity is not rel ated to any system- 
atic periods often present in ASM data (see iFarrell et al.l 120051 : 
Benlloch 2004). We speculate that this period is caused by the 
orbit of the neutron star around its Be companion. It is very close 
to the lower limit of 23 d as proposed by Koyama e t al, (199 lb , 
assuming a circular orbit for the binary system. If the 1997 and 
2002 outbursts are removed from the lightcurve, the 20.85 d pe- 
riod is still present, albeit with a reduced significance. The fact 
that inclusion of the outburst data boosts the signal is expected 
from a Be/X-ray binary, where it is common for outbursts to oc- 
cur near periastron passage of the neutron star However, if the 
orbital period of this system is 20.85 d then it is clear that we do 
not see bright outbursts at every periastron passage. After its ini- 
tial discovery in 1972, Cep X-4 went into outburst every four to 
five years (1988, 1993, 1997, 2002). If this is a continuing trend, 
we may expect to see an outburst in 2007. 

4. Spectral Analysis 

All observations from Table [1] were added to create the average 
outburst spectrum. 

In the 3.5-70 keV energy range we used a power law with 
Fermi-Dirac cutoff to describe the continuum. This model has 
the analytical form 

FDCOiE)^AE-^^^-^±-^ (1) 

and differs from the high energy cutoff model in that the tran- 
sition from power law to exponential is smooth, i.e. without a 
discontinuity of the derivative a t the cutoff en ergy. Although the 
high energy cutoff and NPEX (lMiharalll995]) models also pro- 
duced adequate fits to the continuum, the FDCO model provides 
both a lower xl over the energy range and smaller uncertainties 
on the model parameters. 

The cutoff energy was fix ed at 17keV from the fit with the 
lowest as determined by iKovama et al.l (1199 ll) during the 
1988 outburst of Cep X-4. This approach was deemed reason- 
able, as preliminary fits to the data converged at values close to 
this cutoff energy. In addition we modeled Fe Ka emission with a 
Gaussian emission line, took into account photoelectric absorp- 
tion using the phabs model in XSPEC, with photoelectric cross- 
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Table 2. Spectral parameters for the average outburst spectrum. 
Errors indicate 90% confidence intervals. 



Parameter 



Value 



r 

■Bait 
■Bfold 

Nu 
Ec 

O-c 
Tc 
Efe 

FeEW 

^Gauss 

Gauss EW 

Flux (3.5-10 keV) 

Xl (dof) 



1 44+0.08 

-0.05 

17keV 
10.9!°;^ keV 
2.r0;2x 10-2cm-2 
30.7+ 1 ;5 keV 
3.6:2.« ]jev 

0.7+02 

63r°fe keV 
0.32+55;"^ keV 
73.3 ±0.1 eV 
14.4+°-2 keV 
4.6*11 keV 
1.8+|^;^ergcm-2 s"' 
0.74(43) 



sections from lBalucinska-Church & McCammonI ( 1992h and el- 
ement abundances from lAnders & Grevessd (Il989l) . We added 
a ~5keV wide emission line at iscauss ~14keV to account for 
variation in the continuum over this energy range. This emis- 
sion line is used to account for a feature noticed in spectra of a 
nu mber of accreting X-ray pulsars, e.g. GS 1843+00 or Her X- 
1 (ICoburnll20bll) . It is not an instrumental feature, as it has 
been observed with RXT E, Ginga (Mi hara.1995) and BeppoSAX 
jSantangelo et al. I I1998I) . and should be included in future mod- 
els of pulsar X-ray continua. An absorption line with a Gaussian 
optical depth profile (Coburnetal. 2002, Eq. 6) was used to 
model a weak cyclotron resonant scattering feature (CRSF) at 
30.7 keV. The F-test probabihty of this impr ovement being by 
chance is 1 X 10^"*. See lProtassov et al] ( l2002l) . however, for lim- 
itations of the F-test in these circumstances. The average out- 
burst spectrum is shown in Fig. [3] while the model parameters 
are given in Table |2l 

The average outburst spectrum is strongly influenced by the 
longest (20 ks) observation on MJD 52450.7 . In order to see how 
well this spectrum represented individual observations we nor- 
malised the average outburst model to each individual observa- 
tion by multiplying it by a constant. This allows us to notice 
broad changes across the spectrum. In earlier data (e.g., obser- 
vations 01-00, 01-01) the model strongly overestimates the data 
at the soft end (< 7 keV) of the spectrum. Whereas for later ob- 
servations (e.g. 02-02, 02-03) the soft end is underestimated (see 
Fig- El. To investigate this effect we fitted the spectra of individ- 
ual observations, this time allowing the normalisation parame- 
ters, the power law index and photoelectric absorption parame- 
ter to vary. Our fits show a clear steepening (from 1.19 + 0.06 to 
1 .50 + 0.04) of the power law as the outburst progresses, but the 
absorption column, although poorly constrained by the model, is 
consistently high (~ 2.5 ± 0.3 x 10^^ cm^) through the outburst. 
This suggests that the spectrum becomes softer as the outburst 
progresses. There is, however, little variation around the energy 
range containing the cyclotron feature. 

As cyclotron features are heavily dependent on the viewing 
angle and geometry of the line-forming region, a phase resolved 
study of Cep X-4 would be beneficial to understanding this sys- 
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Fig. 4. Residuals, plotted as Ax, as a function of energy for the 
average outburst model compared to the spectra of observations 
01-00 and 02-03. 




Fig. 5. Pulse profiles as a function of energy for observation 01- 
00 (MJD 52439.9) shown in counts/s/PCU. Profiles are plotted 
twice for clarity. Hardness ratio is 10-20keV band divided by 
the 2-5 ke V band. 

tem. However, the low count statistics in the region around the 
cyclotron line in this dataset make such a study unfeasible. 

5. Pulse Profiles 



An epoch folding search dLeahv et al.ll 19831) was used to deter- 
mine the source pulse period for each observation (see Table[T]i. 
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Fig. 3. The spectrum of Cep X-4 from 3.5~ 
70 keV. The crosses show the data, the smooth 
curve shows the unfolded spectrum, and the 
histogram the model fit. In the second panel the 
residuals, plotted as A^, are shown for the case 
in which no CRSF is included in the model. In 
the lower panel, these residuals are shown again 
- this time including a CRSF in the model. 



The pulse profiles were generated by folding the lightcurves at a 
period of 66.30 s. The pulse profiles are double-peaked and com- 
plex through all energy bands (Fig. |5]l. Both of the peaks show 
additional structure which is prominent at low energies and be- 
comes less coherent at higher energies. The pulse profile in the 
2-5 keV energy changes significantly as the source luminosity 
drops (see Fig. |6]l with the second peak becoming weaker rela- 
tive to the first as the outburst progresses. 

The pulsed fraction, as defined by (F^ax - ^min)/^max, is 
shown in Fig. [T] and indicates an overall increase in softer 
bands (2-5 keV and 5-lOkeV) as the outburst progresses. This 
is contrary to what was observed in the 1997 outburst, where 
IWilson et al.1 (11999 ') saw the pulsed fraction decreasing as the 
source luminosity dropped. 



6. Discussion 

6.1. Spectrum 

A cyclotron line is identified in the average outburst spectrum 
at an energy of 30.7^[ gkeV, confirming the cyclotron line dis- 
covered at t h is ene rgy in the 1988 outburst of Cep X-4 by 
iMihara et aP (Il99lh with Ginga. In the 2002 outburst, as in 
1988, the cyclotron feature appears as an absorption line. With 
the intention of comparing the line properties between the 2002 
and 1998 outbursts we fi t the 2002 data with the CYCLABS 
model, which was used bv lMihara et al.l (Il99lb . The line param- 
eters for the 2002 outburst, using the CYCLABS model are con- 
sistent, within errors, with those measured using the Gaussian 
optical depth model and given in Table |2] However, the overall 
model fit is slightly worse. We find values of the line depth and 
width that are around a factor three smaller than th ose measured 
in 19 88 (Wc = 15.0+ 1.4 keV and De ^ 2.93±0.10. lMihara et al.1 
11991 ). but note that the cyclotron line was near the upper limit 
of the Ginga energy range. The observed energy of the cyclotron 
line is given by 



Assuming that the observed feature is the fundamental (n - I) 
and assuming a gravitational redshift z = 0.3 at the surface of a 
typical neutron star mass of 1.4 and radius of 10 km, we can 
obtain the neutron star magnetic field to be B ^ (3.4 + 0.2) x 
lO'^Gfrom Eq|2] 

The Fe Ka line has an equivalent width of 73 eV which is 
within the range of the Fe line equivalent width during the 1988 
outburst, and is of the same order as that of Fe Ka emission lines 
observed in other accreting X-ray pulsars (Nagase 1989)- 

The A^H column density is more than double that noted be- 
fore in any previous outburst of Cep X-4 (Kovama et al.l ll99Tl : 
Schulz et al. 1995) and also exceeds that derived fro m the red- 
dening of the optical counterpart feonnet-Bidaud & Mouche3 
ll998i) . We interpret this increase in the column density as a local 
effect and ascribe it to possible increased mass transfer through 
the accretion stream or warping of the accretion disc. 

The spectrum certainly changes with source luminosity, be- 
coming softer as the luminosit y fades. A similar effect oc- 
curred in the 1988 outburst with Kovama e t al.l ( 11991 ') reporting 
a change in the power law index from 1.10 + 0.01 to 1.14 + 0.01 
over 10 days of the fading outburst. Although the change in 
power law slope occurs over a similar timescale in this latest 
data set, the change is more pronounced. This is in contrast to 
V0332-I-53, an accreting X-ray pulsar where the spectrum be- 
came harder through the decline of the outburst ( Mowlavi et aT| 
2006). More observations, including data below 3keV, will be 
useful in assessing whether both the power law and H i absorp- 
tion change over the outburst. 

6.2. Pulse Profiles 



£c ^ 11.6 keVx 



1 



B 



l+z 1012G 



(2) 



For the 1997 outburst Wils on et al.l (19991) describe the pulse 
profile using the model of B rainerd & Meszarosl(ll99lh in which 
the total pulse profile is made up of two components. A hard 
pencil beam is caused by photons from one magnetic pole be- 
ing upscattered through the accretion column. These photons 
contribute to the prominent peak in the hardness ratio around 
phase 0.3. Simultaneously, a soft, double-lobed peak originates 
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Fig. 6. Pulse profiles in the 2-5 keV energy for a number of ob- 
servations throughout the outburst. Luminosity decreases from 
the top to bottom panel. (See Table 1 for luminosity values.) 



from the antipodal magnetic pole when cyclotron photons are 
backscattered. These photons are backscattered into the neutron 
star or gravitationally focused around the neutron star to form a 
fan beam of soft photons. 

The pulse profiles of Cep X-4 bear a s triking resemblance 
to the complex pulse profiles of Vela X- 1 dKrevkenbohm et al] 
120021) in the soft X-rays. Both these accreting systems show 
two main pulses each subsequently made up of further struc- 
ture. As with Vela X-1 this complexity disappears at higher ener- 
gies leaving a clearly double-peaked profile. Interpretation of the 
complexity is uncertain, rangin g from variable absorption over 
the neutron star spin phase (iNagase et alJll983i) to anisotropic 
emission through the accretion column at the polar cap. 

The pulse profiles, especially those at lower energies, show 
clear evolution over the outburst. Over the 12 days covered by 
the RXTE pointed observations, the luminosity decreased by a 
factor of 5 and the two peaks in the double-pulsed profiles in 
the 2-5 keV band became distinctly more dissimilar to each 
other with the flux in one peak decreasing far more sharply with 
luminosity than the decrease in flux in the other pulse peak. 
Comparable changes w ere noted in the 1997 outburst of Cep X-4 
dMukerjee et aLluOOOh . where the relative strengths of the two 
pulses comprising the double-pulsed profile were reversed and 
the interpulse became stronger with decreasing source luminos- 
ity. 
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Fig. 7. Pulsed fraction through the outburst. Triangles represent 
the pulsed fraction in the 10-20keV band, diamonds for the 5- 
10 ke V band and circles for the 2-5 ke V band. 



It is not uncommon for binary pulsars to show strong varia- 
tions in their pulse profiles over the duration of an outburst. For 
EXO 2030+375, which also shows complex double- peaked pro- 
files, a reversal of the dominant pulse was detected (Parm ar et al.l 
1989). A simple geometric model comprising pencil and fan 
beams from two offset magnetic poles was employed. This 
interprets the above reversal as a switching of the dominant 
radiation from one magnetic pole to the other. iParmar et al.l 
(Il989h also showed that as the luminosity decreased by a fac- 
tor of 100, the beam patt ern changed from fan to pencil beam. 
iBasko & SunvaevI d 19761) propose that a fan beam will be the 
dominant pulse shape at higher luminosities (> lO^'^ergs"') 
where the accreting matter forms a shock above the neutron star 
surface and radiation escapes predominan tly from the side of 
the accretion column (I Wang & FrankI 19811) . In lower luminosity 
scenarios, the infalling matter may be decelerated by Coulomb 
interactions at the neutron star surface (iBasko & Sunvaevll 19751; 
Kirk & Gallowav 1981), giving rise to a pencil beam emission 
pattern. 

The pulse profi les in Fig. [5] and Fig. [6| show behaviour sim- 
ilar to V0332+53 dTsygankovet al.l 120061) in its high luminos- 
ity (~ lO-'^ergs"') states, i.e., a change in the relative heights 
of the double peaked profile with increasing energy as well as 
a marked change in the relative intensity of the peaks in the 
soft bands (3-6 keV and 6-10 keV) as the source luminosity de- 
creases. Luminosities > 10^^ erg s"' are certainly not seen in this 
outburst of Cep X-4, yet it shows similar changes of the pulse 
profile. This may indicate that the soft pulse profile evolution 
in this source could be attributed to changes in the flux con- 
tribution from the two magnetic poles, although the prominent 
feature at phase 0.3 in the hardness profile (Fig. |5]) is not obvi- 
ously explained by such a scenario. Accurate orbital parameters 
for Cep X-4, which are not yet determined, together with model- 
ing of the pulse profiles will help to highlight the geometry and 
clarify our understanding of the pulse profile behaviour 



7. Summary 

1. We have confirmed the detecti on of a cyclotron line first 
noted by iMihara et aTl dl99ll) in the 1988 outburst of 
Cep X-4. 
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2. We have observed a column density a factor of two higher 
than previously observed for this source. We attribute this 
increased density to an effect local to the Be/X-ray binary 
system of Cep X-4, such as possible partial obscuration by 
the accretion stream onto the neutron star. 

3. We note changes in pulse profiles both with energy and with 
decreasing source luminosity. Although similar changes are 
noticed in the pulse profiles of other accreting X-ray pul- 
sars there exists no global interpretation of these effects. 
Modeling of individual systems, incorporating the orbital pa- 
rameters, can shed light on the pulsar geometry and emission 
patterns. 

4. A softening of the source spectrum with decreasing lu- 
minosity, as was no ted in a previous outburst of Cep X-4 
jKovamaet allliggi t) is observed, in contrast to the harden- 
ing noticed during the outburst decay of V0332H-53. 

5. A tentative orbital period of 20.85 d for Cep X-4 is revealed 
in the long term X-ray lightcurve. 

Acknowledgements. VAM would like to acknowledge the NRF (S.Africa), the 
British Council and Southampton University. 
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